Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 29 April 2010 (MN MgX style file v2.2) 



Evolution of isolated neutron stars till accretion. 
The role of initial magnetic field 



P.A. Boldin \ S.B. Popov 2 * 

1 National Research Nuclear University " MEPhI" , Kashirskoe shosse 31, Moscow, 115409, Russia 
2 Sternberg Astronomical Institute, Universitctski pr. 13, Moscow, 119991, Russia 



Accepted Received ; in original form 



ABSTRACT 

We study evolution of isolated neutron stars on long time scale and 
calculate distribution of these sources in the main evolutionary 
stages: Ejector, Propeller, Accretor, and Georotator. We compare 
different initial magnetic field distributions taking into account a 
possibility of magnetic field decay, and include in our calculations 
the stage of subsonic Propeller. 

ft is shown that though the subsonic propeller stage can 
be relatively long, initially highly magnetized neutron stars 
(Bo > 10 13 G) reach the accretion regime within the Galactic life- 
time if their kick velocities are not too large. The fact that in previ- 
ous studies made >10 years ago, such objects were not considered 
results in a slight increase of the Accretor fraction in comparison 
with earlier conclusions. Most of the neutron stars similar to the 
Magnificent seven are expected to become accreting from the in- 
terstellar medium after few billion years of their evolution. They 
are the main predecestors of accreting isolated neutron stars. 

Key words: stars: neutron — pulsars: general 



1 INTRODUCTION 

Accreting isolated neutron stars (AINS ) were 
predicted 40 years ago by Shvartsmanl J1971I - ) 
and independently by lOstriker^T^iT (|l970h . 
In early 90s there was some enthusiasm due 
to the launch of the ROSAT satellite, which 
was e xpected to fin d man y sources of this 
kind (jTreves fc Cofpil Il99lh . Seve ral popula- 
tions ! studies have bee n ma de (|Blaes et al 

19901 ; iBlaes fc Raiagopall [l99ll: iBlaes fc Madau 



jop 

es 



19951: IMadau fc Blaed 119941 ; IBlaes et all Il995f : 
Manning et al.lll996t ). However, it came out that 
AINS, if they exist, are very elusive (Colpi 
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et al. 1998). The main reason is that initial 
(kick) velocities of NSs appeared to be signif- 
icantl y larger, than it hav e been thought be- 
fore jLvne fc Lorimei] [l994). Initial guess that 
the number of Accretors is small due to low lu- 
minosity of high- velocity NS s was shown to be 
wrong. In a detailed study bv lPopov et"ail (2000) 
(hereafter Paper I) it was shown that INS with 
Crab-like initial parameters and constant mag- 
netic fields spend all their lives as Eje ctors (we 
follow the classification summarized in iLipunovl 
1992) if their initial velocities are > 100 km s _ . 
Then, the fraction of Accretors was mainly deter- 
mined by the fraction of low-velocity NSs. 

Up to the very end of 90s, it was believed 
that the wast majority of NSs are born sim- 
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ilar to the Crab pulsar. I.e., that they have 
short initial spin periods (from milliseconds to 
few tens of millisecond) and magnetic fields 
B ~ I0 12 G. Now it is believed, that about 
one half of NSs have different initial proper - 
ties (|Pqpov et alJl200d : iKeane fc Kramer! l2008h . 
There are at least three groups of sources 
with distinct parameters: compact central objects 
(CCOs) in supernova remnants (SNR), magne- 
tars (anomalous X-ray pulsars - AXPs, and soft 
gamma-ray repeaters - SGRs), and cooling radio- 
quiet NSs du bbed the Magnificent seven (M7) 
(|Popovl |200Sl). CCOs have low initial fields ~ 
10 11 G (jflalpern et al.ll2007l ; iGotthelf fc Halpernl 
2009) and relatively long spin periods (hun- 
dreds of millisecond). AXPs and SGR s have large 
fields ~ 10 14 G (see a review in iMereghettil 
120081 ). The Magnificent se ven-like NSs h a ve fields 
slight ly above 10 13 G (|Haberll 120071 ; iKaplanl 
|2008|) . Probabl y, some of rotating ra dio tran- 
sients fRRATs. McLaughlin et alJlsjQOd ) are sim- 
ilar to the M7. This variety in initial proper- 
ties deserves new studies of evolving NSs using 
th e population synthesis tech nique (see a review 
in iPopov fc Prokhorovl [2007). In this paper we 
present the first step. 

We describe two models. At first, we discuss 
a simple semianalytical approach, which is used 
to illustrate the main features of the scenario. In 
this model velocities and ambient densities are 
not changing. Then we present a detailed numer- 
ical model, which takes into account spatial move- 
ments of NSs in the Galactic potential and real- 
istic 3D distribution of the interstellar medium 
(ISM). Our main results are based on this model. 

In the next section we present basic concepts 
used in both models, and describe each of them. 
Then, in Sec. 3, we present results. Discussion is 
given in Sec. 4. In the last section we present our 
conclusions. 



2 MODELS 

In this section we describe our models. We start 
with explanation of some basic processes and pa- 
rameters of magneto-rotational evolution used in 
both models. Then we discuss the semianalytical 
and the full numerical model, consequently. 



2.1 Basic processes and parameters 

Here we describe some aspects of magneto- 
rotational evolution implemented in both semi- 
analytical and numerical models. 



2.1.1 Standard magneto-rotational evolution 

Here we mainly follow the approach described in 
iLipunovl l|l992h . We consider a NS being born as 
an Ejector. At this stage a relativistic wind and 
Poynting flux are so strong that they prevent in- 
coming matter to penetrate inside neither gravi- 
tational capture radius, Rg, nor inside the light 
cylinder radius, R\. Rg represent the typical scale 
at which the ISM is captured by the NS gravity: 



R G = 2GM/vl 



(1) 



where M is a NS mass and v xe \ is a relative ve- 
locity of the NS and the ISM. The light cylinder 
radius is defined as: 



R\ = c/cj, 



(2) 



where u is the spin frequency of a NS. The ejected 
matter creates a cavern in the ISM within the dis- 
tance of the Shvartsman radius, R s h, at which the 
magneto-dipole pressure, P ~ /j 2 / Rf Rg h , is equal 
to the ram pressure of the ISM, P ~ pv 2 cl . At this 
stage a young NS can be visible as a radiopulsar 
(PSR), and we assume that it losses energy via 
relativistic wind and Poynting flux according to 
the magneto-dipole formula: 



1 dluj 2 

2 dt 



2 /J, id .2 



(3) 



Here I = 10 45 g cm 2 is a moment of inertia of a 
NS, /i is a magnetic dipole moment, % ls an angle 
between rotational and magnetic axis, which is 
assumed to be 7r/2 everywhere below. 

The Ejector stage finish when the Shvarts- 
man radius, P s h, becomes less than Rg- The 
regime changes because after matter appear in- 
side Rg, its pressure start to grow Pmattor ~ 
r -5/2 rQjjg more ra pidiy than the relativistic 

wind pressure growth: P w i n d ~ r~ 2 . So, the con- 
dition Pmattor > Pvind is reached and the pulsar 
wind cannot stop the incoming flow. 

Another reason that causes the Ejector stage 
to cease is disappearence of the magneto-dipole 
emission. This happens when P s h becomes less 
than the light cylinder radius. So, now matter 
fills the light cylinder, preventing the creation of 
the magneto-dipole emission. 
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For isolated NSs both cases can happen not 
only due to spin-down, but also because an ob- 
ject enters a more dense region of the ISM. Vari- 
ations of the relative velocity of a NS and the 
ISM are important, too. For small velocities _Rg 
can become larger than i? s h. Oppositely, for large 
velocities i? s h can become smaller than the light 
cylinder radius due to the ram pressure (in this 
case Rg < Ri due to a large velocity). 

To summarize, when at least one of condi- 
tions i? a h < Rg or i? s h < R\ is met, the Ejector 
stage ceases. After that, matter falls down till its 
pressure is counterbalanced by the magnetic field 
pressure. The radius at which these pressures are 
equalized is called the Alfven radius: 



Ra 



2/7 



for the case Ra < Rg and 

1/6 

Ra 



(4) 



(5) 



for the opposite. M is the accretion rate. 

The Propeller stage begins if Ra > Rc, where 



Rc 



(GM\ 



1/3 



(6) 



is the corotation radius at which the solid body 
rotational velocity equals the escape velocity. At 
this stage matter is "propelled" away from a 
star, because of its interaction with the magneto- 
sphere. 

For the prope l ler sta ge we use the model pro- 
posed bv lShakurall|l975f ). The same approach was 
also used in Paper I. The period derivative can be 
written as: 



dP 
dt 



= MRlPr 1 ~ KP° 



For a = 1 we obtain: 

K = 2.4 x 10- 14 ^V/Vo 9/T s"\ 



(7) 



(8) 



where /i3o is a magnetic dipole moment in units 
of 10 30 G cm 3 , n is the ISM number density, 
Wio is the total velocity \/ a'§ + v 2 cl in units of 10 
km s~\ a B — 10 km s _1 is the sound speed. 

When a NS spin-downs enough for the con- 
dition R c > Ra to be met, a star leaves the 
Propeller stage and switches to the next stage, 
depending on the relation between Ra and Rg. 
Note, that in our scenario this can also happen 
because of changes in the ISM density or in the 
velocity of a NS. 



If -Ra > Rg, then a NS enters the Georota- 
tor stage, which is called so because of similarity 
of the NS magnetosphere structure to the Earth 
magnetosphere in the fast solar wind. At this 
stage in our model no spin-up/spin-down mecha- 
nisms are taken into account. 

Otherwise, if Ra < Rg then a star at 
first enters the subsonic Propeller stage, where 
the main accretion mechanism - the Rayleigh- 
Taylor instability - is supressed by high tem- 
perature of an envelope (i.e., the gas is "too 
light" ) . Temperature increases because of heating 
(a NS loses angular momentum to the envelope 
and heats it), and decreases because of radiative 
losses (mainly bremsstrahlung) . The spin-down 
rate at this stage is tak en in the following form 
(|Davies fc Pringldll98lh : 



dP 
~dt 



2.4 x lO^Vao m" 1 s s" 1 . 



(0) 



Note, that P is constant if the field is not decay- 
ing. 

A NS starts to accrete from the ISM and en- 
ters the Accretor stage when heating of an enve- 
lope due to rotational energy losses by a NS be- 
comes less effective than bremsstrahlung cooling, 
and so the matter at the magnetospheric bound- 
ary becomes so heavy that the Rayleigh- Taylor 
instabil ity develops. T his occurs at the critical 
period l|lkhsanoj20o"if ): 



8.7 x 10 4 R A % M3 2/3 m 1/6 s 



(10) 



where _Ra, io is the Alfven radius in units of 
10 cm. At the Accretor stage the spin-down 
is calculated as at the subsonic Propeller stage 
(eq. [9|, and we neglect possible quasi equilib- 
rium at this stage. Th is is an oversimplification 
(|Prokhorov et al.ll2002T ). but as here we are not 
interested in details of spin properties of Accre- 
tors we neglect some details. 



2.1.2 Magnetic field distribution and decay. 
Initial period distribution 

We consider two different field distributions. The 
first is a delta-function. As a standard value we 
use /X30 = 1, which gives us B cq — 10 12 G, accord- 
ing to fi = B Bq R 3 . This distribution is chosen to 
make comparison with the results from Paper I . 

The second distribution is a result of the 
magnetic field dec ay starting with the "optimal" 
one from the paper lPopov et all (j201Ch . This "op- 
timal" distribution is the lognormal one with 
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<log(B pole /[G])) = 13.25 and cr logSpolc = 0.6, 
where B po i c is the value of the poloidal field on 
the magnetic pole (-B po ic = 2B cq ). Till the mag- 
netic field reaches some saturation value B m m it 
undergoes decay according to 



B(t) = B 



-t/-r oil „ 



T Hall v 



-*/TOJuo) 



(11) 



where rohm is the Ohmic characteristic time, 
and THaii is the typical timescale of the fast, ini- 
tial Hall stage, that depends on the initial field 
(THaii oc 1 /B ). Typically, r hm = 10 6 yrs and 
THaii = 10 4 yrs (for B = 10 15 G). The asymp- 
totic value of the field depends of the initial 
strength. In or der to approxima te the results of 
simulations by I Pons et all (2009) we choose 



ZJ,,,;,, = :nw{ ^,2x 10 13 G 



(12) 



Since, according to eq. (|16[) min(tE) ~ 10 7 yr 
3> Tohm = 10 6 yr, we neglect the process of decay 
and take the field already decayed down to P m i n 
as the inital field in our model. The shape of this 
distribution is shown in Fig. 

Below we will refer to these two distributions 
as "the standard" (for fiso = 1) and "the de- 
cayed" respectively. Everywhere below B is the 
polar field. 

For Po in the complete model we take the 
distribution with (Po) = 0.2 5 s and gp n = 0.1 s. 
Such distribution was used in IPopov et al1(|2010T ). 
In our scenario here this is a simplification, as a 
NS can spin-down while its field decays down to 
the minimum value. So, an initial period in our 
model should be different (longer, depending on 
the strength of the initial magneti c field) from a 
period used in IPopov et all l|2010h . This results 
in slight overestimating of the number of objects 
at the Ejector stage on the price of other stages. 
But we tested that this assumption does not influ- 
ence our results significantly even for magnetar- 
scale initial fields. This is so because a star always 
evolves off the Ejector stage much slower than the 
field decays. 



defined for any evolutionary stage, and it just 
demonstrates how efficiently a NS interacts with 
the surrounding medium. Almost all characteris- 
tic radii and critical periods depend on M. 

Not only the absolute value of a kick is im- 
portant, its direction is significant, too. If a kick 
has large component perpendicular to the Galac- 
tic disc, w z , then a star spends much less time 
close to the Galactic plane, where the ISM den- 
sity is higher (see below). 

Most of Accretors in our scenario have small 
velocities, so only low-velocity end of the distri- 
bution is important. Several shapes of the kick ve- 
locity di stribution have been discu ss ed in the lit- 
erature: Arzoumanian et al. (l2002l ): )Hobbs et all 
l|2005h : iFaucher-Gieuere fc Kaspi (|2006l l Here 
we use the initial velocity distri bution pro- 
posed by lArzoumanian et all l|2002h . It is a bi- 
maxwellian distribution with the Gaussian three- 
dimensional dispersions <7i = 90 km s _1 and 
oi = 500 km s" 1 . We vary the contribution of 
each of these components using the parameter w\ , 
which is the fraction of NSs in the low-velocity 
component. 



2.2 Semianalytical model 

In a simple semianalytical model we assume 
monotonic magneto-rotational evolution of a NS 
in constant conditions: NS velocity, its magnetic 
field and the ISM density do not change during a 
NS lifetime. Each transition from stage to stage is 
defined by solving appropriate equations for the 
given characterisic values. In this model we al- 
ways use the NS mass M = 1.4Mq (m = 1.4) 
and the initial period Po = 0.02 s. 

At first, we solve equations P s h(P, fl,v) = 
Rg, P s h(P, v) = Pi(P) to find the period at 
which the Ejector stage ends and the star be- 
comes a Propeller. This gives us: 



P(E-^P)«7^« 
for R B h = Pg and 



1/2 1/2 -1/4 



-1/2 



(13) 



2.1.3 Velocity distribution 

A NS initial velocity (and so, its kick) has great 
impact on its magneto-rotational evolution. This 
is due to the fact that the efficient accretion rate 
M strongly depends on the velocity: M ~ v^Jf. 
Note, that M (i.e. just a combination of the ISM 
density, relative velocity and NS mass) can be 



P(E ±% P) 



i An 1/3 -1/3 -1/6 

142 fi 3 ' v 10 n s 



(14) 



for P s h = Pi. The critical period (for any branch 
of the transition from the Ejector stage) is called 
below Pe 

Integrating eq. @ we obtain: 



P(i) = 



3c 3 7 



(15) 
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Then, neglecting Po we derive the time of the first 
transition: 

t(E -A P) « 8.25xl0V3o 1 ^io«" 1/2 ^" 1 yrs, (16) 
t(E P) « 3.26xl0 1 V3o 4/ Vo 2/E n" 1/3 yrs.(17) 

Among these values the smaller one is used. 

If for a given fi and v the value of £e exceeds 
teal — 10 10 yr, then the fraction of lifetime which 
a star spends as an Ejector, te, is equal to 1. 
Elsewhere, te is equal to tE/tGai- 

Next, we must consider the Propeller stage. 
Solving the equation R C (P) = Ra(^,v) for both 
- Ra > Rg and Pa < Rg - cases, we obtain the 
critical periods: 



P(P A A G i»P) = 500A$'i$'r» 



6/7 9/7 -3/7 -11/7 /, Q \ 
' m ' S, (18) 



P(P 



A>G 



G) = 3xl0 5 At^ 2 «r 1/2 «~ 1/4 ^~ 1/2 s.(19) 



These periods correspond to the end of the Pro- 
peller stage. Any of such periods is called below 
Pp (note, that before in several papers such a pe- 
riod was called Pa, as without the subsonic Pro- 
peller stage and neglecting the possibility that a 
star becomes a Georotator, after reaching Pp ac- 
cretion starts). 

Then by solving the equation log (Pa) — 
log(PE) = KAt on At we obtain Q the value for 
Atp. This is the time period during which a NS 
stays at the Propeller stage. The fraction of Pro- 
pellers in the total distribution among evolution- 
ary stages is taken in the form rp = min(l — 
te, Atp/iGai). 

If i?G(u) < Ra{h,v) < Pc(Pp) a star is con- 
sidered as a Georotator. The fraction of NSs at 
this stage is tg = 1 — (te + rp). 

If Ra{/J-,v) < Pg(«), and due to the spin- 
down at the Propeller stage Ra((J,,v) < P c (Pa), 
then the subsonic Propeller stage begins. The 
spin-down during this stage is given by eq. (|9j). 
The stage ends when P = Pb re ak, which gives 
us the equation: At = (Pb rC ak — Pa) /P, where 
P is taken according to eq. (|9j), from which we 
obtain Ai ss p - the duration of the subsonic Pro- 
peller stage. Fraction of this stage is r ss p = 
min(l - (te + rp) , At ssP /t G ai) 

After a NS spin-downs and leaves the sub- 
sonic Propeller stage, it starts to accrete. The 
fraction of Accretors is ta = 1 — (te + rp + r ss p). 



1 Pe is taken as P(E — y P) for the case Ra < Pd 
for the opposite case. 



and P(E i% P) 



For the scenario without the subsonic Pro- 
peller stage we assumed r ss p = (or, equivalenlty, 
taking ta = 1 — (te + rp) , as it was done in Paper 
I)- 



2.3 Complete numerical model 

Using the same subroutines for characteristic val- 
ues and spin-down rates as in the semianalytical 
model, we make a more detailed numerical model. 
A NS evolution in this model proceeds in realis- 
tic conditions: the Galactic potential and the ISM 
density distribution. 

This is done by splitting all the time from 
to tcai into an equdistant grid. Then, on the 
domain of the acquired time grid we compute all 
values which do not independ on the rotation of 
a star: velocities, coordinates, ISM densities, Pa, 
M, and so on. 

After all rotation-independent values are 
fetched, the magneto-rotational evolution is cal- 
culated starting with the Ejector stage. Transi- 
tions between stages are determined by changes 
in relations between characteristic values, as de- 
scribed above. Note, that there is almost no pro- 
hibited transitions, because of large variations in 
the environment and in the velocity of a NS. 



2.3.1 Spatial evolution 

Magneto-rotational evolution strongly depends 
on the spatial evolution of a NS, moving through 
the Galaxy. A NS recieves an initial kick velocity. 
The vectors of the kick and the progenitor's Ke- 
plerian velocities are summed. In our model, the 
ISM rotates with a Keplerian velocity around the 
center of the Galaxy. 

Equations of motion of a NS with a given 
initial values are solved on the time grid us - 
ing the LSODA subroutine (|Hindmarshl Il983h . 
The Galactic potential is taken in the same form 
as in Paper I (i.e. in the form suggested in 
iMivamoto fc Nagailll975l and Paczynski 1990). It 
is a three component potential (disc, buldge, and 
halo) which reproduces well enough trajectories 
on a long time scale (billions of years). In the 
problem concidered here we do not need a more 
complicated gravitational potential. 
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2.3.2 ISM density and NS initial spatial 
distribution 

The ISM density is taken according t o the old an- 
alytical model in lPosselt et al.l (|2008l ). It is gener- 
ally the same as in Paper I with some corrections 
in the z-dependence (see below). This distribution 
have exponential or Guassian behaviour perpen- 
dicular to the Galactic disc. The radial distribu- 
tion has a peak at R ~ 5 kpc. For very low-density 
regions (large R and z) we used the minimum 
value of the ISM number denstity n = 10~ 5 cm -3 . 

In Paper I there w as a small mistake (copied 
from lZane et al.lll995h regarding the ISM distri- 
bution. Dispersions in eqs. (5) and (7) of Paper 
I should be divided by 2.35 (as they are actually 
not dispersions, but FWHM). In eq. (6) of Paper 
I the coefficients 0.345, 0.107, and 0.064 should 
be 0.7, 0.19, and 0.11, correspondently. 

Exactly as in Paper I, the birthrate of NSs is 
proportional to the square of the local ISM den- 
sity. 



3 RESULTS 

3.1 Semianalytical model results 

With this model we address two main questions. 
How does the account for the subsonic Propeller 
stage influence the fraction of Accretors? How 
does the fraction of accretors depend on the mag- 
netic field? 

In Paper I the subsonic Propeller stage was 
not used, and only Crab-like fields, ~ 10 12 G, 
were considered. Here we include this stage and 
take into account higher magnetic fields (up to 
the values typical for decayed fields of magne- 
tars). Obviously, the first effect reduces the num- 
ber of Accretors, while the second - increases. It 
is interesting to understand with a simple model 
the interplay between them before addressing the 
same questions with a more advanced one. In this 
model we use only the delta-function magnetic 
field distribution. 

In Fig. 1 we show the dependence of the frac- 
tion of INSs at the Accretor stage on the magnetic 
field and the ISM number density. Here both of 
these parameters are constant during a NS life- 
time. For Crab-like fields and the ISM density 
~ 0.1 - 1 cm" 3 it is below ~1%. The num- 
ber of Accretors for a given realistic ISM den- 
sity steadily grows with increasing initial mag- 
netic fields up to ~ 5 x 10 13 G. This is related to 




Figure 1. Fraction of Accretors in the semianalyti- 
cal model. The bi-maxwellian kick velocity distribu- 
tion from Arzoumanian, Chcrnoff and Cordes (2002) 
is taken. In the left panel we show results for the sce- 
nario with the subsonic Propeller stage. In the right - 
we neglect this stage to demonstrate its influence on 
our results. 

the fact that a NS with larger magnetic field spin- 
downs faster, and so quicker reaches the stage of 
accretion. However, for larger fields it happens at 
large periods, but the increase of the spin-down 
rate is more important than the increase of criti- 
cal periods for transitions. 

For large fields the fraction of Accretors 
starts to decrease, because more and more sources 
appear as Georotators due to their huge magneto- 
spheres. However, it is important to repeat, that 
here we assume constant fields, but it is normally 
accepted that large fields gradualy decay while a 
NS is aging. 

The difference with Paper I, where for /j, « 
10 30 G cm 3 , n fa 1 cm -3 and comparable ve- 
locities we obtained ~ several percents of Accre- 
tors, is explained by the influence of the subsonic 
Propeller stage. It is visible in the right panel. 
The fraction of accretors for n = 1 cm -3 and 
fi — 10 30 G cm 3 is about one percent if the sub- 
sonic Propeller stage is not taken into account. 
For the initial field in the M7 range the fraction 
of Accretors in our new model (with the subsonic 
Propeller stage included, see the left panel) goes 
up to ~ 10 - 30%. 

Clearly, the main conclusion here is that the 
subsonic Propeller is not a strong barrier for NSs 
with reali stically large field s in contrast with con- 
clusion by Ikhsanov ( 2001]) . As now we know that 
the fraction of such objects (SGRs, AXPs, M7, 
RRATs) is not low - tens of percent, - one can ex- 
pect that significant fraction of old INS can start 
to accrete. Thus, one has to study the distribu- 
tion of INSs in different evolutionary stages in 
more details. 

Some results of this subsection are summa- 
rized in the Table. 1. We demonstrate twelve 
tracks for different fi, n, v. NSs following tracks II, 
III, and IX always stay at the Ejector stage. High 
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Table 1. NS evolutionary tracks 



Track 


n, cm a 


M30 


via 


TE 


Pe, s 


TP 


Pp, 


3 


TssP 


P break ) S 


Track I 


0.5 


1 


5 


0.419 


16.051 


0.423 


3.163 x 


10 3 


0.850 


2.278 x 10 6 


Track II 


0.5 


1 


20 
















Track III 


0.5 


1 


40 
















Track IV 


0.5 


10 


5 


0.042 


50.758 


0.042 


2.276 x 


10 4 


0.067 


1.317 x 10 7 


Track V 


0.5 


10 


20 


0.168 


101.517 


0.170 


1.353 x 


10 6 


0.651 


2.568 x 10 8 


Track VI 


0.5 


10 


40 


0.163 


100.091 


0.169 


1.523 x 


10 5 


Georotator 


Track VII 


2.0 


1 


5 


0.209 


11.350 


0.212 


1.746 x 


10 3 


0.370 


8.464 x 10 5 


Track VIII 


2.0 


1 


20 


0.838 


22.700 


0.854 


1.038 x 


10 4 






Track IX 


2.0 


1 


40 
















Track X 


2.0 


10 


5 


0.021 


35.892 


0.021 


1.257 x 


10 4 


0.030 


4.892 x 10 6 


Track XI 


2.0 


10 


20 


0.084 


71.783 


0.085 


7.469 x 


10 4 


0.264 


9.541 x 10 7 


Track XII 


2.0 


10 


40 


0.103 


79.442 


0.106 


1.077 x 


10 5 


Georotator 



velocity NSs with fields larger than the Crab- 
like value, become Georotators after the Propeller 
stage (tracks VI and XII). The NS following track 
VIII never becomes an Accretor as it stays for a 
long time at the subsonic Propeller stage. All the 
rest NSs (tracks I, IV, V, VII, X, XI) finally start 
to accrete. 



3.2 Complete numerical model 

In this subsection we present our main results 
obtained with the complete model for 10 5 evolu- 
tionary tracks. Our main aim is to calculate the 
distribution in evolutionary stages for two differ- 
ent distributions of the initial magnetic field. In 
addition, we demonstrate the effect of changing 
velocity distribution. 

As we use the bi-Maxwellian distribu- 
tion proposed by lArzoumanian et al.l (|2002h . to 
demonstrate the dependence of our results on 
the velocity distribution, we decided to change 
relative contributions of the two constituents. In 
FigEH he horizontal axis shows w\ - the contribu- 
tion of the low- velocity part of the bi-Maxwellian 
distribution. For w\ = we have a pure 
Maxwellian distribution with a = 500 km s _1 , for 
Wi = 1 — a pure Maxwellian with a — 90 km s - . 

We show results of calculations for two distri- 
butions of initial magnetic fields described above 
(sec. 2. 1.2). The first is just a delta-function n = 
10 30 G cm 3 . It corresponds to the typical assump- 
tion made in 90s. The second is based on recent 
results bv lPopov et al.l (|2010| ). 

Fractions of Ejectors, Propellers, subsonic 
Propellers, Accretors and Georotators demon- 
strate monotonic, nearly linear behavior. The 
number of Ejectors strongly decreases with in- 



creasing wi. The behavior of Accretors and sub- 
sonic Propellers is opposite. 

The behavior for the two studied field distri- 
butions is similar in the cases of Accretors, Ejec- 
tors, and subsonic Propellers. In the case of Pro- 
pellers and Georotators the situation is different 
for two distributions. 

In the most realist ic case according to 
lArzoumanian et al.l l|2002h . - wi = 0.4 - we have 
(in the case of initially decayed field distribu- 
tion) ~ 55% of Ejectors, ~ 5% of supersonic and 
~ 20% of subsonic Propellers, ~ 10% of Accre- 
tors, and finally, ~ 10% of Georotators. As we 
see, now for our "the best choice" model we pre- 
dict more Accretors than in Paper I. It is what 
was expected on the basis of the semianalytical 
model. 

The increase in the relative number of Ac- 
cretors is due to the presence of INSs with large 
initial magnetic fields. This is illustrated in Fig[3] 
We show there contributions of INSs with dif- 
ferent initial magnetic fields to the population 
of Accretors. Note, that the scale is logarithmic 
in both axis. INSs with initial fields < 3 10 12 G 
are more numerous than those with 10 13 G < 
B < 2 10 13 G0 However, the latter produce seven 
times more Accretors. Still, many (about 1/2) of 
NSs with the largest initial field considered here 
do not produce many Accretors as they become 
Georotators due to large spatial velocities. 

Typically, NSs become Accretors in the re- 
gions of the Galactic disc, where the ISM den- 
sity is higher. As it was noted before NSs with 



2 Here we speak about decayed fields in our model. If 
we recalculate it to obtain real pre-decay initial fields, 
then the intervals are changed. 
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Figure 2. Different panels show fractions of INSs 
at different stages in the complete numerical model: 
Accretors (top left), Propellers (top center), sub- 
sonic Propellers (top right), Ejectors (bottom left) 
and Georotators (bottom center). Black circles de- 
note the decayed magnetic field distribution, whereas 
empty squares corresponds to the delta-function fi = 
10 30 G cm 3 . Bottom right panel shows the ratio of 
Propellers to Georotators. Triangles correspond to 
an additional variant of magnetic field distribution. 
In this variant only low-field - B < 5 x 10 12 G - 
INS are taken from the decayed model. By w± we 
mark the fraction of low-velocity component of the 
bi-maxwcllian distribution from Arzoumanian, Cher- 
noff and Cordes (2002). Here fj, = BR 3 /2, as B is the 
magnetic field value at the pole. 

low total velocity but significant z-component, v z , 
spend most of its lifetime outside the Galactic 
plane. Most likely the longest stage for such a NS 
is the subsonic Propeller. If a star is born rela- 
tively far from the Galactic center and recieves a 
large kick then it escapes from the Galaxy. Such 
a NS spends most of its life as an Ejector. Alter- 
natively, if the velocity is high but not enough to 
escape the Galaxy, the NS returns to the Galac- 
tic plane after some long time, then it can quickly 
pass the Propeller stage (so-called non- gravitating 
Propeller in this case) and become a Georotator. 

In addition to the global distribution we com- 
pute separately distributions over stages inside 
(R < 16 kpc and \z\ < 1 kpc) and outside the 
Galaxy (in the following paragraphs we refer as 
"the Galaxy" only to the former volume). It can 
seem suprising and confusing, but according to 
our model we predict more Accretors than Ejec- 
tors inside the Galaxy: ta ~ 30% and te ~ 20%. 
Subsonic Propellers are more abundant than Ac- 
cretors and Ejectors in this volume: r as p ~ 43%. 
These numbers can be explained in the follow- 
ing way. Most of NSs which contribute a lot to 
the number of Accretors have v z < 100 km s _1 . 
They spend most of their lives inside the Galaxy, 
and so there they dominate. Other stages, corre- 



spondently, are not abundunt inside the Galaxy: 
Georotators contribute tg ~ 7%, Propellers - 
Tp ~ 1 — 2%. Roughly, NSs with kick veloci- 
ties from the low-velocity part of the distribution 
(about 30-40 %) stay inside the Galaxy. Those 
with magnetic fields higher than typical radio pul- 
sar values become Accretors. 

The situation outside the cylinder R < 16 
kpc and \z\ < 1 kpc is the following. Ejectors 
contribute te = 76%, Georotators - tg ~ 11%, 
subsonic Propellers - r as p = 9%, Propellers - 
Tp w 4%, and Accretors - ta ~ 0.2%. As one 
can see, the situation with Ejectors and Accre- 
tors is opposite in comparison with the internal 
part. Almost all Accretors are situated inside the 
Galaxy. Note, that due to this Fig. [3] with the 
total distribution also refers to the population of 
Accretors inside the Galaxy. Almost 2/3 of all 
neutron stars are outside the Galaxy and they 
are either Ejectors or Georotators. 

In the solar neighborhood (i? so iar < 2 kpc 
and \z\ < 0.5 kpc) we predict ~ 35 - 40% of Ac- 
cretors and slightly more (~ 40 — 45%) subsonic 
Propellers with only ~ 18 — 20% of Ejectors. Con- 
tributions of others stages are negligible. In to- 
tal, in the solar proximity (-R so i ar < 2 kpc and 
\z\ < 0.5 kpc) there are 0.33% of all NSs. This 
gives us, for A*ns = 10 9 the number density in 
the solar neighborhood no ~ 3 x 10~ 4 pc ~ 3 , in 
good correspondence with recent results bv lOfekl 
(2009) and with earlier studies. 



4 DISCUSSION 

Available estimates of the number of Accretors 
should be taken with care, since there are sev- 
eral effects which act in inhibiting accretion, as it 
is discussed below. I.e., the number of observable 
Accretors can be much smaller. Still, low velocity 
INSs moving through a high density medium (the 
only ones with high enough luminosities to be 
potentially detected) and with strong magnetic 
fields can become accretors after < few Gyr. It 
was shown above that INSs like the M7, which 
have magnetic fields higher that those typical of 
radio pulsars, and at the same time have not very 
large spatial velocities, are the most favored as 
Accretors predecessors. 

Although the number of AINSs might be 
even larger than that originally estimated in Pa- 
per I, the conclusion that at low fluxes Accretors 
outnumber cooling isolated NS (Coolers) is based 
on the assumption that the luminosity corre- 
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Figure 3. The dashed line represents the field dis- 
tribution after the immediate field decay discussed in 
sec. 2.1.2, i.e. the decayed model. The solid line shows 
the contribution to the number of Accretors for dif- 
ferent initial fields in the decayed model, normalized 
to unity. The velocity distribution is calculated for 
uti = 0.4. Here B is the magnetic field value at the 
magnetic pole. Each bin shows mean fraction of time, 
among all the stars having such a fields. 



sponds to the Bon di accretion rate. T his is a quite 
controversial issue. iBlaes et all (|l995t l have shown 
that for typical ISM densities (n S3 1 cm -3 ) ac- 
cretion rate does not exceed ~ few 10 9 gs _1 , even 
if the star velocity drops below ~ 60 kms -1 . This 
is due to the ionization of the ISM surrounding 
the star by the X-ray radiation which, in turn, 
produces an increase in the sound speed freez- 
ing the accretion rate. However, in Paper I it 
was shown that the velocity distribution of Ac- 
cretors peaks at ~ 50 kins -1 , and for these ve- 
locities the effect is small. So, for most Accretors, 
heating of the ISM can be neglected, especially 
if they appear in regions of high ISM density. A 
further issue is the role played by the star mag- 
netic field in the accretion flow dynamics outside 
the Prope ller stage. On th e basis of 2D MHD cal- 
culations, iToropina et"afl l|2003h concluded that 
only a fraction of the initial (Bondi) flow reaches 
the star surface, and this fraction decreases with 
growing magnetic field of a NS. Whether 3D in- 
stabilities may counteract this effect is still an 
open question. 

Anyway, if several weak sources without 
measured proper motions and interpreted as INS 
candidates can be identified in ROS AT, Chandra 
or/an d XMM- Newton archives fsee lTurner et all 
l20ld and references therein), or discovered by 



eROSITA, then it is not trivial to distinguish 
Coolers from Accretors. 

If an INS comes to the Accretor stage 
only after a long subsonic Propeller episode, 
its spin period is > 10 4 s. Such long periods 
are not unexpected even if the subsonic Pro- 
peller stage is neglected. When a NS starts 
to accrete it continues to spin down, until it 
reaches a quasi-equili brium period, P ea sa 10 6 s 
for n = 1 cm" 3 (jKonenkov fc Popovl 1 19971 : 
IProkhorov et al.ll2002r i. The ultra-long spin peri- 
ods of Accretors could be the best discriminator 
between this type of sources and Coolers, which 
are expected to have spin periods < few seconds 
(like the M7 and cooling PSRs). However, at low 
fluxes it would be extremely difficult to discover 
pulsations in Coolers, so the non-detection of a 
periodicity is not a strong argument in favour of 
an AINS. 

Opposite to Coolers, Accretors are expected 
to show both - spin-up and spin-down - as their 
periods fluctuate around the quasi-equilibrium 
value. However, P measurements can be impossi- 
ble for faint sources with very long periods. 

The period of accreting INSs can be signifi- 
cantly shorter than ~ 10 6 s in the case of mag- 
netic fields decaying down to small values (~ 
10 9 G), although some kind of fine tuning is nec- 
essary. As discussed above, to reach accretion in a 
time shorter than the Hubble time an INS should 
have at least a magnetic field S3 10 12 G. So, decay 
should not be significant during the first ~ 1 Gyr 
of the evolution, otherwise a N S spends all its life 
as an Ejector or a Propeller ( C^l^LetaL 19981 : 
iLivio et alJ Il998l : IPopqv fc Prokhorovl l2000h . If 



the field decays during the Accretor (or even sub- 
sonic Propeller) phase, an INS can attain a period 
~ 10 3 -10 4 s, since P cq is smaller for smaller fields. 

Accretors, at variance with coolers, are not 
expected to be steady sources because of changes 
in the accretion rate, due to inhomogeneities of 
the ISM, on a time-scale 



/ _ Rg „ 1n 8 -3 

t ss ~ 3 x 10 v 10 s . 

V 



(20) 



Note that this time scale is shorter for fainter 
sources. 

Spatial distribution of Accretors and new 
weaker Coolers are expected to be slightly dif- 
ferent, as the first represent much older popula- 
tion, and for the first higher ISM density is fa- 
vorable for detection in contrast with the second. 
N ew (i.e., undiscovere d, yet) Coolers according 
to lPosselt et alJ lj200ct ) are expected to be found 
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at distances ~ 1 kpc. So, they should be rela- 
tively bright, ~ 10 31 erg s _1 . Accretors cannot 
be that bright, and so they are expected to be 
found closer. Young Coolers should trace star- 
forming regions. Accretors, which already expe- 
rienced long evolution in the Galactic potential, 
should be distributed more smoothly. However, 
for them to be detectable it is important to be 
inside regions of relatively high ISM density. 

The X-ray spectrum of a NS accreting at low 
rate from the ISM is very similar to those of cool- 
ing INSs, at lea st in the case when the latter has a 
H atmosphere ll Treves et al.ll200ol and references 
therein). Nevertheless, for the same luminosity, 
the effective temperature of an Accretor is higher 
and, hence, the spectrum is harder because of sig- 
nificantly reduced emitting area. For typical val- 
ues of the star mass and radius, the hot polar cap 
size is 



D n c w ln3 -2/7 -3/7 1/7 

-Reap ~ 9.5 x 10 fi 30 v 10 ' n 



(21) 



This is smaller than the size of a typical emitting 
area in Coolers. Spectra of Accretors are expected 
to be harder than those of Coolers. 

We summarize some differences between Ac- 
cretors and cooling NSs (Coolers) in Table 2. 



5 CONCLUSIONS 

After the first of the M7 have been discov- 
ered (|Walter et al.l Il99rj ). several authors pro- 
osed and di s cusse d that they can be AINSs 



posed and di scussed t hat they can be A1IN as 
(I Walter et all Il99fj; iKonenkov fc Popovl Il997l ; 
iNeuhauser fc Triimpei 19991 ). Though, it ap- 
peared that it is not so. The M7 are young NSs 
with relatively large fields. P robably, they are re - 
lated to evolved magnetars (|Pqpov et al.l 120101 ). 
Here we demonstrate that in future the M7 and 
similar sources are expected to become AINS if 
their magnetic fields do not decay significantly. 
Even a relatively lon g stage of subsonic Pro- 
peller (|lkhsanovll200lh cannot prevent accretion. 
This is a good news for observers. Probably, tele- 
scopes like eROSITA aboard Spektr-RG will be 
able to detect AINS, soon. However, the ques- 
tion of the accretion efficiency is still on the list 
(|Toropina et alJl2003h . 

The distribution over evolutionary stages 
strongly depends on kick velocity distribution, 
initial magnetic field distribution and field evolu- 
tion. Because of that precise predictions are not 
possible now. This shows how important is to de- 
tect old isolated NSs as Accretors (or, less proba- 



ble, other stages) to learn more about initial prop- 
erties and evotuion of INSs. 
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Table 2. Comparison of properties of dim Accretors and Coolers 





Accretors 


Coolers 


Spectrum 


Harder, ~ hundreds eV 


Softer, ~ 100 eV 


Spin periods 


Very long, > 10 s s 


Shorter, ~ 10 s 


V 


Variable 


Stable spin-down 


Distance 


Close, ~ 100 - 200 pc 


Further away, ~ 1 kpc 


Luminosity 


Low, ~ 10 29 erg s — 1 


Higher, ~ 10 31 erg s" 1 


Variability 


Variable, At ~ weeks - yrs 


Stable 


Spatial distribution 


Towards higher gas density 


Towards starforming regions 
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